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ABSTRACT 

Blue straggler stars (BSS) are abundantly observed in all Galactic globular clusters (GGC) where data exist. 
However, observations alone cannot reveal the relative importance of various formation channels or the typical 
formation times for this well studied population of non-standard stars. Using a state of the art Henon-type 
Monte Carlo code that includes all relevant physical processes, and observationally motivated initial conditions 
in large ranges in density, concentration, binary fraction, and mass, we create a little over 100 models with 
properties typical of the observed GGCs. These models can be directly compared with the observed GGCs 
and BSSs can be easily identified. We identify the BSSs in our models and determine their formation channels 
and times of creation. We find that for central densities above 10 3 M Q pc~ 3 the dominant formation channel 
is stellar collisions. The majority of these collisions happen via binary-mediated strong scattering encounters. 
As a result a strong correlation between the total number of BSSs and the binary fraction in a cluster can be 
observed. We find that the number of BSSs in the core show weak correlation with the collision rate V. Close 
inspection of our models reveals that the observed lack of correlation between the number of BSSs and V for the 
observed GGCs may be because the estimate of V using observed cluster properties is not a reliable indicator 
of the true T, Our models indicate that the BSSs observed today may have formed many Gyrs ago. Denser 
clusters have higher contributions from collisions which typically create higher mass BSSs compared to those 
created via mass-transfer in a binary. As a result, we find that the BSSs in denser clusters typically have formed 
more recently than those in sparser clusters. 

Subject headings: methods: numerical — methods: statistical — blue stragglers — stars: kinematics and 
dynamics — globular clusters: general 



1. INTRODUCTION 

Blue straggler stars (BSSs) are stars bluer and brighter (and 
therefo re more massiv e) than the main-sequence tumoff (in a 
cluster: [Sandagd 1 95 3l) . It is impossible for normal single star 
evolution to create them. Instead, blue stragglers must be cre- 
ated when a main-sequence (MS) star is disturbed either via 
mass transfer in a binary or by a collision with another star. 
Both processes can increase the stellar mass and possibly also 
rejuvenate the s tar by bringing in new Hydrogen for burning 
in the core (eg.,[Lom bardi et ali i 19951 1 19961: ISills et al.lll997l 
2001; Lombardi et al. 2002h. Since the BSSs are among the 
brighter members of a cluster, they are relatively easier to de- 
tect in photometric observations. All Galactic globular clus- 
ters (GGCs) where the data exi st have shown this exotic stellar 
popul ation in abundance (e.g., Pio tto et alJ20 02: Ferra roet alJ 
1995). Many open clusters als o show moderately larg e popu- 
lations (~ 10s) of BSSs (e.g. JMathieu & GeUeil206^) . 

A large body of theoretical work has been devoted to ex- 
plaining the formation and properties of BSSs. This work 
suggest that there are two different formation channels for 
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the BSSs, namely col lisions involving at least one MS star 
(Lom bardi et al.lfl995h and mass transfer in a binary (MTB) 
from a more evolved star (giant or another MS star ) to a MS 
star (for a detailed recent study IChen & Hani 120091) . In both 
cases the key physical effect is replenishing the MS star with 
fresh Hydrogen (via a collision or via MTB), and increasing 
the star's mass. If this Hydrogen can be mixed inside the MS 
star and be brought to the core of the MS star for burning, then 
the MS star lives in the Hydrogen burning stage longer than 
an undisturbed star of equal stellar mass is expected to live. 
The star is then called rejuvenated. This rejuvenated star then 
would be observed as a BSS when its mass is larger than the 
standard turnoff mass in the cluster. 

Stellar collisions can happen in various ways in a cluster. 
Two single stars can physically collide (SS collision). If bina- 
ries are present in the cluster, then binaries interact with other 
single stars (BS) or binaries (BB). During the BS and BB in- 
teractions collisions can oc cur between two or more stars tak- 
ing part in the interactions (Fregeau et al. 20QH). 

In mass transferring binaries the details of the orbital evolu- 
tion and mass transfer rate determines its final outcome. Mass 
transfer can be stable, where, a fraction of the companion 
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mass is transferred to the MS star. In contras t, dynamical 
mass transfer leads to merger of the two stars dChen & Hani 
2009). Stellar mergers may also happen in interactions in- 
volving higher order hierarchical systems. For example, a hi- 
erarchical triple under certain configurations can incre ase the 
eccen tricity of the inner binary orbit via Kozai effect dKozail 
|1962|) . If the eccentricity is sufficiently hi gh, the inner binary 
stars m ay be merged through tidal friction (Perets & Fa brvckvl 
2009). Although it may be a signi ficant channel of BSS for - 
mation in less dense open clusters (iPerets & Fabrvckvll2009l) . 
it is not clear whether these triples can be stable in more dense 
regions typical of the GGCs and thus be as effective. 

Since BSSs are MS stars that have been disturbed, and pos- 
sibly rejuvenated (by some channel), some form of interac- 
tion must have created them. These interactions are in turn 
dependent on the cluster global properties, including the cen- 
tral density (p c ), binary fraction (/;,), and mass-segregation 
timescale. Thus the BSSs in a cluster can in principle pro- 
vide a window into the cluster's dynamical history. In this 
respect the formation channel(s) of the BSSs and their typical 
formation times are of great interest. Unfortunately, although 
observations in various wavelength bands have taught us a lot 
about the stellar properties of the BSSs, it is hard to determine 
how a particular BSS formed or how long ago the interaction 
occurred. Theoretical modeling is the only way to learn this. 

A popular approach so far to learn about the observed BSSs 
in the GGCs has been by carefully studying the correlations 
between the number of BSSs (normalized in some way by an- 
other characteristic stellar population in these clusters) with 
some dynamically importa nt cluster property, i ncluding the 
collision ra te estimator T (|Poplev & Hutl 120061) . the binary 
fraction ft (Sollima et al. 2008), and total cluster mass (most 
recently ILeigh et aLll2013l) . While these observed correlations 
can teach us a lot about the possible formation channels of 
these exotic stellar populations, direct understanding of these 
processes are impossible to achieve. For example, a BSS cre- 
ated via mass transfer in a binary could actually be a result 
of some past dynamical encounter that had changed the or- 
bital properties of a binary to alter its normal evolutionary 
path in such a way that mass transfer occurs. On the other 
hand, a collision in fact may be mediated by resonant dy- 
namical interactions involving stellar binaries. In addition to 
these complicated interactions, it is also hard to know obser- 
vationally exactly what the binary fraction in the cluster is 
now and is impossible to know what it had been in the past. 
Similarly, the calculation of the T parameter involves many 
properties in the cluster that are hard to determine. The only 
way to explore and probe these complicated interactions in an 
evolving cluster is by actual A-body modeling of a realistic 
cluster with realistic A, stellar mass function, and binary frac- 
tion and including all relevant physical processes such as two- 
body relaxation, strong encounters including physical colli- 
sions, binary-binary (BB) and binary-single (BS) interactions, 
stellar evolution, and Galactic tidal stripping. In these mod- 
els if all encounters are documented in addition to the over- 
all cluster properties and stellar properties at all times, then 
all necessary information can be obtained. BSSs can then be 
identified from a snapshot at any given time from these mod- 
els and for them the full history from t = can simply be 
looked up. 

Until recently realistic modeling of GGCs was hard to do 
in a star-by-star detail including all relevant physical pro- 
cesses. Hence, idealized simulations have been used to study 
the formation channels and their radial distribution in the clus- 



ter dMapelli et al.ll2.Q06h . However, Northwestern group's re- 
cently modified Henon-type Monte Carlo code CMC is well 
suited to quickly model GGCs including all relevant physi- 
cal processes. Using this code it is now possible to model 
many clusters with star-by-star detail using a large range of 
observationally motivated initial parameters. For example, a 
recent study with CMC has explored a large range in param- 
eter space and successfully created a population of ~ 200 re- 
alistic cluster mo dels that have properties ver y similar to the 
observed GGCs dChatteriee et al.ll2010i [20T3b . These mod- 
els include all relevant physical processes, and since they also 
have realistic A, these models can be directly compared with 
the observed GGCs without any need for scaling. In addition 
to the dynamical properties, single and binary stellar evolu- 
tion properties are also followed in these models. Thus these 
simulations provide an unprecedented opportunity to under- 
stand the dominant dynamical processes responsible for BSS 
production in typical GGCs. In addition, since these models 
also span a large range in dynamically important global prop- 
erties typically observed in the GGCs, these models also can 
shed light on the widely used correlations between the num- 
ber of BSSs in a cluster and the cluster's global properties and 
show what these correlations truly mean. 

This study is the first in a series where these models are 
explored, BSSs are extracted, their dynamical history is ob- 
tained and two very basic questions are answered. What is 
the dominant physical process that creates the BSSs observed 
in the GGCs? What is the typical age since formation for the 
BSSs? 

This study is organized as follows. In Section |2] we de- 
scribe briefly our simulation setup, describe how the BSSs are 
identified in the models, and definitions for various categories 
of BSSs based on their formation channels. In Section [3] we 
show examples of the radial distribution for our model BSSs. 
We present the relative importance of the various formation 
channels in Section H The time elapsed since formation for 
the BSSs in our models is explored in Section [5] In Section 
[6] we investigate correlations between the specific frequency 
of the BSSs in the core with /& and V. We summarize and 
conclude in Section [7] 

2. NUMERICAL METHOD 

We use our Henon-type Cluster Monte Carlo code (CMC) 
to numerically model star clusters typical of the observed 
GGCs including realistic number (A) of stars and primordial 
binary fraction (ft,). We include all relevant physical pro- 
cesses namely two-body relaxation, strong interactions such 
as physical collisions and binary mediated scattering interac- 
tions, Galactic tidal stripping, and single and binary stellar 
evolution. This code is developed and tested rigorously over 
the past decade (iJoshi et alJl2000i feOOlUFregeau et aljbool 
i Fregeau &~ Rasio 2007; Chatte riee et al.l 12.010b lUmbreit et all 
l2012t iPattabiraman et al.ll2013t IChatteriee et al.ll2013l) . Us^ 
ing a large grid of initial conditions over the range of values 
typical of the observed young massive clusters we create 121 
detailed star-by-star models and evolve them for ^12 Gyr. 
Our choice of initial conditions is based on observations of su- 
per star clusters. All our simulated clusters have initial virial 
radii between r v =3-4 pc. Initial A is varied between 4 
- 10 x 10 5 stars. The positions and velocities of the stars 
(and center of mass for the binaries) are assigned according 
to King profiles with the dimensionless concentration param- 
eter Wq between 4-8. We vary the initial binary fraction, ft, 
between 0.05 — 0.3. The stellar masses of the stars, or pri- 
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m aries in case of a binary, are chosen from the IMF presented 
in lKro upa (2001, their Equations 1 and 2) in the stellar mass 
range 0.1 — 100 M Q . Secondary binary companion masses 
are sampled from a uniform distribution of mass ratios in the 
range 0.1 M© — m p , where m p is the mass of the primary. 
The semi-major axes, a, for stellar binaries are chosen from 
a distribution flat in log a within physical limits, namely be- 
tween 5 x the physical contact of the components and the local 
hard-soft boundary. Although initially all binaries in our mod- 
els are hard at their respective positions, some of these hard 
binaries can become soft during the evolution of the cluster. 
The cluster contracts under two-body relaxation and the ve- 
locity dispersion increases making initially hard binaries soft. 
Moreover, binaries sink to the core due to mass segregation 
where the velocity dispersion is higher than that at the initial 
binary positions. We include these soft binaries in our simu- 
lations until they are naturally disrupted via strong encounters 
in the cluster. 

The numerical setup and initial properties for these simu- 
lation s are similar to the work presented in (Cha tterjee et al.1 
2013) which show that these models at the end of the simula- 
tion (at cluster age t c \ w 12 Gyr) attain properties including 
the core radius (r c ), the half-mass radius (r/j), central density 
(p c ), total mass (M c \), and relaxation timescale at half-mass 
(t r h) very similar to those observed in the GGCs. The only 
difference between t hese simulations and those presented in 
IChatterjee et al.1 (1201 3l) is that models with zero initial /(, are 
discarded in this study since observation s suggest that this 
is not the case fo r the GGCs in reality dDavis et al.l 12008b 
iMilone et alJl2012l) . In addition, the models presented in this 
study explores a larger range in initial /&, and a slightly larger 
range in initial King concentration parameter Wq . Using these 
highly realistic models we focus on understanding and char- 
acterizing the BSSs in globular clusters. For a detailed com- 
parison of model properties wit h the observed proper ties of 
the GGCs readers are referred to IChatterjee et al.l(l2013t . 

2.1. Definition of blue straggler stars in our models 

We define BSSs as stars that are still in the Hydrogen burn- 
ing stage but with masses > 1.1 A/to where, Mro is the 
main-sequence turn-off mass of the cluster. Note that in real 
clusters BSSs are chosen from a HRD as stars that are brighter 
and hotter than the MS tum-off. Unfortunately, the defining 
box for the BSSs vary from cluster to cluster depending on, 
for example, the quality of the data, and the scatter width of 
the MS (e.g., Leigh et al. 2007). Often the criterion used 
by the observers is that no MS star should be counted as a 
BSS and to ascertain that some actual BSSs can be discarded 
because of their proximity to the MS in the HRD. Thus it is 
hard to formulate a self-consistent luminosity and temperature 
based criterion that would be valid for all clusters. The mass- 
based criterion is much simpler and self-consistent to use for 
the models (although, inaccessible to observers). We find that 
in our models the BSSs identified using the mass-based crite- 
rion populate the expected region of synthetic HRDs created 
using our models (Figure Q] for an example). Since in this 
study we restrict ourselves to the theoretical understanding of 
the BSSs, in particular to understand the dominant formation 
channels and typical ages since formation, we simply adopt 
the mass-based definition for BSSs. 

Due to the differences between the BSS identification cri- 
teria in our models and for observed clusters we find that al- 
though our models populate similar ranges in global proper- 
ties such as r c with the observed GGCs, we can identify more 



BSSs than observed. For example, Figure|2]shows the number 
of BSSs within the core as a function of the core radius both 
for our models and for observed GGCs. The observed BSS 
counts are taken from the compilation in iLeigh et al.l (120071) . 
Our models clearly have very similar core radii compared to 
the observed GGCs, however, we identify a larger number 
of BSSs in the core. This is expected since due to the un- 
certainty in determining the main-sequence turn-off in real 
clusters a significant number of BSSs can remain unidenti- 
fied due to their proximity to the MS. Indeed, as part of a 
future study we have seen that if a more observationally moti- 
vated selection crite ria for our model B SSs is used following 
the prescriptions of Leig h et all (120071) . the number of BSSs 
match well with the observations. In addition, the numbers 
of BSSs using either criteria are well correlated (Sills et al. 
2013, in preparation). In addition to the differences in the 
identification criteria for the BSSs in our models with the ob- 
served, one should also note that within BSE, the rejuvena- 
tion prescription used for mergers, collisions and binary mass 
transfers ass umes "full mixing" (for details see iHurley et al.l 
I2000ll2002l) . This is a simplification. In reality, the amount 
of Hydrogen mixing and thus the degree of rejuvenation of 
the M S star depends on the k inematic details of the encounter 
(e.g., ISills et al.|[l997l |2001l) . Because of this simplification 
the lifetimes of the BSSs after rejuvenation in our models are 
actually upper limits for non-rotating stars. This also leads 
to an over-production of BSSs in our models to some degree. 
This issue is further discussed in Section [5] 

2.2. Classification of blue-straggler stars based on their 
formation channels 

At the last snapshot of the model (t c \) we identify the 
BSSs following the criterion described in Section |2] The 
full dynamical history for each of these BSSs is retrieved. 
We identify the interaction (such as SS collision, BS colli- 
sion, or MTB) that had increased the mass of this star for 
the first time such that it would be identified as a BSS at 
t c \. We classify each BSS into one of 5 groups accord- 
ing to this first interaction. The 5 groups are: stable mass 
transfer in a binary denoted as MTB, binary-stellar-evolution- 
driven merger/disruption denoted as SE, single-single colli- 
sions denoted as SS, collisions mediated by binary-single in- 
teractions denoted as BS, and collisions mediated by binary- 
binary encounters denoted as BB. Note that some BSSs may 
have a complicated history. For example, after a mass trans- 
fer episode in a binary, a MS star may have attained a BSS 
mass; later on this star could still collide with another MS 
star, further growing in mass and getting rejuvenated for a 
second time. For BSSs that are formed via such multiple 
channels, only the first one is considered for the classification 
to avoid ambiguity. In this particular example the formation 
channel for this BSS would be classified as MTB. Neverthe- 
less, such complicated multiple episodes of rejuvenation of 
the same MS star are not frequent. There is another source 
of complication. A BSS classified as MTB (created via MTB) 
may have had interacted dynamically prior to the Roche-lobe 
overflow that had not changed its mass. Such an interaction 
would change the orbital properties of the binary, sending it 
along a very different evolutionary pathway that led to the 
mass transfer event that ultimately created the BSS. These 
perturbed MTB or SE channels are duly noted and discussed. 
However, interactions of this type do not change the formation 
classification of the BSSs. 
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FIG. 1 . — Example of a synthetic Hertzsprung-Russell diagram for a representative model with many blue straggler stars (runl 10, Table [T}. Dots denote 
all stars in the cluster. Red squares and blue triangles denote blue-straggler stars in single and binary systems, respectively at the integration stopping time of 
t c l = 12 Gyr. Blue straggler stars are defined as H-burning stars with stellar masses > 1.1 M^o, where Mto is the main-sequence turn-off mass. BSSs in 
this model populate the expected region in the synthetic HRD. In addition, the single and binary main-sequences, giant branch, white-dwarf cooling sequence are 
clearly visible. 



3. RADIAL DISTRIBUTION OF THE BLUE STRAGGLER STARS 

One of the observationally accessible properties of the 
BSSs in the GGCs is their radial distribution. The radial dis- 
tribution of the BSSs in a cluster depends on a collection of 
the host cluster's properties including p c , relaxation time, and 
fb. Thus analysis of the radial distribution of the BSSs has 
been proposed as a tool to observationally constrain the dom- 
inant formation channel, and various dynamical properties of 
the host GGCs. As a result, the radial distribution of the BSSs 
has been st udied with interest for many GGCs via observa- 
tions (e.g.. iFusi Pecci et all 1 19921: iDalessan dro et al.l [2 008b; 
Ferrar oet al.l Il995t iBeccari et al.l 120081: ID alessandr o et al.l 
2008at iMoretti et all l2008allbh IFerraro et al.l 1201 2b as well 



via theoretical modeling (e.g., [Mapelli et al. 2006). 
Three qualitatively different types of radial distributions 
are observed among the GGCs, namely, single-peaked, 
double-peaked, an d fiat distributions (for a compilation see 
IFerraro etal1l20 121) . 

To find the relative frequency of the BSSs in the cluster 
as a function of their radial positions we defi ne a double- 
norma lized BSS frequency i?ess following IFerraro et al.l 
(2012?) given by 



Rbss 



SNbss/N: 



BSS 



SL/L C 



(1) 



where SNbss denotes the number of BSSs within some radial 
bin, TVbss denotes the total number of BSSs, 5L denotes the 
stellar luminosity within that bin, and L c \ denotes the total 
luminosity of the cluster. i?sss is a measure of whether the 
BSSs follow the luminosity density profile of the cluster or 
they are over/under abundant in some parts. 

Figure [3] shows -Rbss as a function of the radial position 
in units of the core radius. As in some observed GGCs (e.g., 
M55) we find examples of double-peaked distributions (bot- 
tom panel). We also find examples of single- peaked dis- 
tribut ions (top panel; e.g., observed in M80, IFerraro et al.l 
120121 and references therein). In our models we do not find 
flat radial distributions obs erved, for example, in NGC 2419 
( IDalessandro et al.l l2008bb . This is expected since all our 
models are relaxed and hence a non-segregated r adial distribu- 
tion fo r the BSSs in our models is not expected. Ferra ro et al.l 
d2012l) suggest that the minimum distance r m j n from the cen- 
ter where the central peak of i?BSS en ds in units of r c is 
anti-correlated with the core relaxation timescale due to mass- 
segregation effects of the relatively high-mass BSSs. Results 
from our models are consistent with this finding. However, we 
find that the determination of r m j n can have large errors and 
depends strongly on the radial bins used to calculate -Rbss- 

4. RELATIVE IMPORTANCE OF BLUE STRAGGLER FORMATION 
CHANNELS 
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FIG. 2. — Number of BSSs in the core (_/Vbss c ) vs r e . Black circles and red squares denote model and observed BSSs, respectively. The observed BSS 
numbers in the core are obtained from Leigh et al. (2007). Our models populate very similar r c values compared to the observed GGCs. However, our models 
show higher Nggg c numbers due to the mass-based selection criteria. In an observed cluster many of these BSSs will not be identified as BSSs due to scatter 
in the main-sequence and the MS turn-off. 



In this section we investigate the relative contributions from 
various BSS formation channels. Figure |4]shows an example 
of the branching ratio for the BSSs in the whole cluster and 
in the core for model run32 (TableQ]). This is a model with 
a moderate central density p c = 7 x 10 3 M Q pc~ 3 . Colli- 
sions produce most of the BSSs in this model. Among the 
collisional BSSs channel SS contributes the least. The bi- 
nary mediated collisional channels, namely BS and BB con- 
tribute about 36% and 50%, respectively. BSSs produced via 
the channel MTB contributes about 11% and SE contributes 
only 3.2%. As expected collisional channels dominate even 
more in this cluster if only core BSSs are considered. Another 
interesting aspect is that in this particular example model al- 
though the overall contribution from collisions via BS and BB 
channels increase if only core BSSs are considered compared 
to all BSSs in the cluster, contribution from the SS channel 
decreases. This trend is seen often in our models. This can 
be counter-intuitive to some degree since due to the higher 
central stellar densities, all forms of stellar collisions are ex- 
pected to be more frequent within the core. This is because of 
the low density of single MS stars within the core at late times 
due to mass segregation. 

About 60% of all BSSs formed via MTB or SE in this exam- 
ple model have had at least one prior strong encounter which 
did not make it a BSS. Thus for this particular model a signifi- 
cant fraction of even the MTB and SE BSSs have been dynam- 



ically modified before binary stellar evolution could rejuve- 
nate them. It will be very interesting to investigate in a future 
study whether such dynamical encounters typically enhance 
the possibility of MTB and SE BSS production relative to the 
field or inhibit it. 

Branching ratios for various formation channels depend di- 
rectly on the global properties, including the central density 
p c of the host cluster. The higher the central density of a clus- 
ter, higher is the chance of a collisional formation of BSSs. 
Figure|5]shows the branching ratios for SS, BS, BB, MTB, and 
SE formation channels as a function of p c . For models with 
central densities p c > few xl0 3 M©pc~ 3 collisional chan- 
nels clearly dominate. In lower density clusters contribution 
from MTB increases. Throughout the full range in binary frac- 
tion in our models (/& = 0.05 - 0.30), the same trend is ob- 
served. This indicates, that even if /& is increased (within the 
studied range in p c ) the relative contribution from MTB does 
not increase proportionally. While increasing ft increases the 
probability that a BSS can be created via MTB (simply be- 
cause there are more binaries), the increased number of bina- 
ries also increases the BS and BB BSS formation rates. Typi- 
cal p c values in the GG Cs are orders of magnitude higher than 
- 10 3 M Q pc- 3 (e.g.. iHarrisI [l99l IChatteriee et al.l I20TI 
I20T31) . Hence, we conclude that the collisional channels dom- 
inate BSS production for the populations observed in major- 
ity of the GGCs. Only about 30% of all GGCs have observed 
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FIG. 3. — Examples of two qualitatively different radial distributions for our model BSSs. Vertical and horizontal axes denote the double-normalized BSS 
frequency Rbss (Equation[T} and radial position in units of r c , respectively. Top panel shows an example of centrally peaked radial distribution (runl 8; Table 
[TJ. Bottom panel shows an example of a double peaked distribution (runl 1 0) one near the center r/r c < 5 and a second peak at large r/r c fa 22. A zone of 
avoidance is observed at r/r c between about 5 and 17. 



all BSSs core BSSs 




FIG. 4. — Example of branching ratios for the various formation channels of BSSs at t = 12 Gyr (run32; Table[T}. This model has a moderate central density 
(p c ~ 7 X 10 3 M0pc — 3 ; Table[T}. This particular model was chosen as an example due to its high number (312) of BSSs. The dominant formation channel 
is stellar collisions in this model. Among the collisional BSS production channels collisions mediated by BB interactions dominate followed by BS mediated 
collisions. SS collisions are rare. Mass transfer in a binary is the dominant channel among the binary-stellar-evolution driven BSS production channels. Binary 
stellar evolution driven mergers or disruptions contribute to a small degree. 
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FIG. 5. — Branching ratios of various formation channels for all BSSs in our models vs the respective cluster's p c . The circles, squares, and triangles denote 
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p c values c onsistent to have BSS populations predominantly 
from MTB (lHarrisll 1 996h . If only the core population of BSSs 
is considered, the importance of the collisional channels in- 
creases as expected. We find that throughout the full range of 
central densities in our models collisional channels dominate 
the formation of core BSSs (Figure [6). 

Although binary stellar evolution is not the dominant chan- 
nel for a large range of central densities typical of the GGCs, 
binaries themselves are nevertheless crucial for significant 
BSS production. Figure |7]shows the relative contribution be- 
tween various collisional formation channels of the BSSs as 
a function of p c . For the full range in /{, in our models BSSs 
produced via SS are rare. Binary-mediated collisions (BS or 
BB) dominate. As f^ increases, so does the relative contribu- 
tion from BB. For example, for models with initial /& = 0.05, 
the BS channel contributes more to the collisional BSSs com- 
pared to the BB channel for all p c values. In contrast, for 
fb = 0.3 the contribution from the BB channel dominates for 
all p c values we have considered. 

Figure [9] shows the fraction of BSSs produced via binary 
stellar evolution (MTB and SE collectively) that have been 
dynamically perturbed before the mass transfer process has 
occurred. These interactions can potentially alter the orbital 
properties of these binaries initiating very different evolution- 
ary pathways that lead to mass transfer. BSSs formed from 
these disturbed binaries clearly would not have formed the 
BSSs with the exact same properties and at the same times 
without these encounters. As expected, for the models with 



higher p c values MTB and SE BSSs have a higher chance of 
being perturbed before production. For a given p c a range of 
these values are also seen. This is due to the different concen- 
tration values. Models with initially lower concentration leads 
to a smaller fraction of perturbed MTB and SE BSSs. The ma- 
jority of the models with p c > 10 4 M Q pc -3 produce more 
perturbed MTB and SE BSSs than unperturbed. Low density 
(p c < 10 3 ) models on the other hand will have a much larger 
fraction of BSSs via the MTB and SE channels created from 
the primordial undisturbed binaries. Note that this is in di- 
rect contrast from so me earlier observational studies, e.g., by 
Solli ma et al.l (12008, further discussion in Section|6}. 

5. TIME SINCE FORMATION FOR THE BLUE STRAGGLER STARS 

How long a BSS will live after the encounter that formed 
and possibly rejuvenated it depends on the supply of new Hy- 
drogen that can be mixed into its core for burning. The de- 
gree of rejuvenatio n depends on the de tailed kinematics of 
the encounter (e . g. . ISills et al J 1997112001b as well as the stel- 
lar properties such a s rotation, and the a ges of the parent stars 
prior to interaction (IChen & H an 2009). Incorporating these 
details is hard and beyond the scope of this study. Instead 
we adopt a full mixing prescription used in the stellar evolu- 
tion code BSE dHurlev et all 120001 12002b that has been used 
in this study as part of our cluster evolution code CMC and 
is widely used in most iV-body codes. It has been previ- 
ously found that this full-mixing prescription adopted in BSE 
predicts a longer rejuvenated lifetime than found in detailed 
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FIG. 7. — Branching ratios between various collisional BSS formation channels vs p c . Squares, triangles, and circles denote BS, BB, and SS collisions, 
respectively. Throughout the studied range of p c and /{, single-single collisions are rare compared to binary-single and binary-binary collision channels. At high 
ft = 0.3 values binary-binary collisions dominate over binary-single collisions for the production of BSSs. 



stellar evolution models assuming head-on collisions (e.g., 
iGlebbeeket ail2 008; Gle bbeek & Polsl2008l) . Hence the time 
since formation for the BSSs in these models (especially cre- 
ated via head-on collisions) should be treated as upper limits. 

We define the elapsed times since formation i agc = t c \ — 
£farm> where tf olm is the time when for the first time a star 
achieved a mass that will give it a BSS status at the final 
cluster age t c \ via some channel. Recent detailed numerical 
studies indicate that most GGCs (observationally categorized 
as the non-core-colla psed clusters) are still contracting (e.g., 
IChatteriee et al.l [20 1 3l) . Hence, relevant global cluster prop- 
erties such as p c , r c , and do not attain steady-state values 
typical of, e.g., core-collapsed a.k.a. binary-burning clusters. 
A priori, a BSS can have any t ag0 between zero and the clus- 
ter's age. However, there may also be a "typical" i ago corre- 
sponding to a peak in the age distribution. Thus to understand 
the environment that had produced the present-day BSSs one 
must take into account the evolving properties of the cluster 
for at least this typical t agc . Unfortunately, t agc is impossible 
to determine purely from observations and thus results from 
detailed cluster models is the only way. 

Figure [10] shows examples of the distributions of t age in 
our models differing by almost two orders of magnitude in 
their p c values. The BSSs in the high-density model are typ- 
ically younger compared to those in the low-density model. 
The broadly peaked distributions seen here are qualitatively 
similar in all our models. However, the exact positions of 



these peaks depend on the cluster properties. For the low- 
density model the bulk of the BSS population have ages be- 
tween about t agc is 0-5 Gyr. This is consistent with pre- 
vious re sults suggested for the observed BSS population in 
47 Tuc (ISills et al.ll2000l) . In contrast, the BSSs in the high- 
density model typically have t agc between 0-3 Gyr. In each 
case there is a long tail in the t agc distribution extending to 
about 6 Gyr. Half of all BSSs in the low-density model is con- 
tained within an age of 2 Gyr, whereas, for the high-density 
model this value is 0.9 Gyr. 

The difference in the peaks of the age distributions between 
the high and low-density models is expected. In a high- 
density model the dominant formation channel for BSSs is 
stellar collisions (SS, BS, and BB; Figures |4][5J. In contrast 
MTB is the dominant formation channel for sufficiently low 
density (p c < 10 3 M0pc -3 ) clusters (Figure |5). Stellar col- 
lisions typically form higher mass BSSs than those formed 
via MTB. Since a lower mass BSS has a longer residual life- 
time than a higher mass counterpart, the typical t agc becomes 
shorter as p c increases due to higher contribution from colli- 
sional BSS-formation channels. Indeed we find that the me- 
dian for the t agc values for the BSSs in our models negatively 
correlate with p c (Figure fTTT i. Due to the same reason the 
median value of t age for the core BSSs is typically slightly 
lower compared to all BSSs. This effect is more prominent 
for higher models. The number of BSSs outside r c is low 
for low fi, models resulting in little difference between the 
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median £ age values for all BSSs in the clusters and for the 
BSSs inside r c only. We find that the oldest BSSs in a cluster 
can be almost as old as the cluster itself. 

6. SPECIFIC NUMBER OF BLUE-STRAGGLER STARS VS BINARY 
FRACTION AND T 

After the large-scale surveys using the Hubble Space 
Telescope (e.g., the ACS survey; ISaraiedini et al.l 120071; 
iDotter et all 120071: iMilone et all 120081) provided high resolu- 
tion observations for the GGCs, a large number of studies 
have investigated various correlations between the number 
of BSSs in these clusters and their gl o bal properties (e.g., 
Leigh et al. 1 120071: ISollima et all 120081 iKnigge et alJ 120091: 
Leigh et al. 2013). In particular, the correlation or lack thereof 
between the core BSSs (A^bss.c) an d two cluster properties 
namely the binary fraction in the core {fb. c ) and the collisional 
T have been studied in much detail. Both of these quantities 
are of great dynamical importance. The presence of primor- 
dial binaries increases the chance of all BSS formation chan- 
nels except SS. Whereas, BSS formation via SS is directly 
related to the collisional T given by 




dPoolev & Hutl [2006). Especially, after it was found that 
the number of low-mass X-ray binarie s and CVs in a clus- 
ter is strongly correlated with T (e.g.. iPoolev & Hulll2006l) 
the same correlation was sought for the BSSs in the core 



since all of them are dynamically created populations. How- 
ever, it was found that the specific number of BSSs in the 
core, normalized by another population similar in mass to ac- 
count for the fact that a higher- N clusters should naturally 
have a higher number of BSSs, e.g., the horizontal branch 
stars, A r BSs,c/A , HB,c, is rather insensitive to the estimated 
r. In contrast, it is found that Nbss,c/Nbb,c shows clear 
correlations with t he core binary fraction in the GGC s (e.g. , 
Leigh et alJ 120071; ISollima et all 120081; IKnigge et alJ 120091; 
Leigh et al]|2013l) . These two pieces of information led to the 
claim that most of the BSSs in the GGCs come via binary stel- 
lar evolution channels (MTB and SE) and collisional channels 
(SS, BS, and BB) are unimportant (e.g., ISollima et"aT1l2008l: 
Knigge et al. 2009). However, our models suggest that colli- 
sions dominate BSS production in majority of the GGCs. We 
now investigate further this apparent discrepancy. 

Figure[l2lshows the specific number of BSSs (A^bss/^Vhb) 
in our models as a function of the overall It also shows 
the specific number of BSSs in the core (A'bss.c/A'hb.c) as a 
function of the core binary fraction fb tC . Clearly, Nbss/Nhb 
is strongly correlated with /& and A^bss.c/^Yhb.c is strongly 
correlated with ft jC . This is not at all surprising. Results in 
Section |4] clearly shows that throughout the studied range in 
p c the dominant BSS production channel is binary-mediated 
collisions (BS and BB). Hence, as fb increases so does the 
number of BSSs. Although, naturally the number of BSSs 
from the MTB and SE channels also increases with increas- 
ing fb, their relative importance is always low compared to 
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FIG. 9. — Fraction of perturbed BSSs created via MTB or SE mergers as a function of p c . Only perturbations before the production of the BSSs are taken into 
account. This fraction is directly related to the fraction of binary stellar evolution driven BSSs (MTB and SE) that would have not been produced with the same 
properties had there not been a previous encounter. In addition, once produced, these BSSs may interact further which was not taken into account. Thus, these 
fractions indicate a lower limit for the fraction of observed BSSs that have not been dynamically altered. Black circles and red squares denote all BSSs in the 
models and only those within the core, respectively. Each panel from top to bottom show initial /j, = 0.05, 0.1, 0.2, and 0.3, respectively. 



BS and BB channels for all p c > 10 3 M pc -3 . Thus the ob- 
served correlation between -/Vbss.c/^Vhb.c with /&. c does not 
actually indicate a predominantly MTB BSS formation, rather 
that most collisions are binary-mediated. 

Figure [T3l shows A^bss,c/-^Vhb,c as a function of T calcu- 
lated using true cluster values in the code. We find a positive 
correlation between the -/Vbss.c/^Vhb.c and the theoretical T 
values. This is only expected since our results suggest that the 
dominant formation channel for the BSSs in typical GGCs is 
collisions of some kind. This appears to be in disagreement 
with the observed trend. Closer inspection reveals that the 
apparent discrepancy is likely coming from the differences in 
theoretically calculated V values and those that are estimated 
from observed clusters. Note that the T values are calculated 
using three cluster properties, r c , p c , and the velocity disper- 
sion (Vc.u) in the core. Most recently Chatte riee et al.l (12013b 
show that it is important to be careful about the definitions of 
all of the above quantities while comparing theoretical models 
with observed. To find what an observer would estimate for 
the r in our model clusters (r o b s to different iate from true T) 
we est imate an observed r c as prescribed in Cha tteriee et al.l 
(120131) . measure the central density as an observer would from 
the surface brightness profile, and calculate v c a taking into 
account only the bright MS (M* > 0.3 Mq) and giant stars in 
the cluster. Using these "observed" values for r c , p c , and v c>(T 
we then calculate the observed value of r o bs for all our mod- 



els. We find that the theoretical T values are not retrieved eas- 
ily from r o bs (Figure \l4\ indicating that the r o b s values esti- 
mated in observed GGCs are not reliable indicators of the ac- 
tual collision rates. Figure [T5lshows ATbss.c/^Vhb.c as a func- 
tion of r o bs for all our models. We find that the correlation be- 
tween ATbss,c/^Vhb,c an d r as seen in Figure [T3lis not appar- 
ent between A/bss,c/^Vhb,c an d r o b s - Thus our models sug- 
gest that although most BSSs are created via binary-mediated 
collisional channels in typical GGCs, A^bss,c/A^hb,c may ap- 
pear insensitive to the observed estimates of r o b s - 

7. CONCLUSIONS 

We have identified the BSSs from a large ensemble of real- 
istic simulated clusters using a simple and well defined mass- 
based criterion. The BSSs populate the expected region in 
a synthetic HRD (Figure [TJ. Although the stellar properties 
and numbers of BSSs in various clusters can be determined 
observationally, it is hard to estimate the typical i ago for ob- 
served BSSs. Moreover, the formation channels of these ob- 
served BSSs are also unknown. Our detailed realistic TV-body 
models provide, for the first time, rigorous constraints on the 
branching ratios between various formation channels, distri- 
butions of t agc and also shed light on the observed correlations 
between the number of BSSs in a cluster and its dynamical 
properties including and T. 

We have identified the relative importance of the various 
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FIG. 10. — Examples for distributions of the time elapsed since formation (t ag e) for the BSSs for two representative models, one for low (run 61) and another 
for high p c (runl 1 3; Table[TJ with central densities varying by two orders of magnitude. The solid (black) and dashed (blue) lines represent the low and high 
density clusters, respectively. The top and bottom panels are for all BSSs and those within r c , respectively. For MTB BSSs the formation time tf orm is set at 
the beginning of Roche-lobe overflow. The BSSs in the low-density model shows a larger age-spread compared to those in the high-density model. Low-density 
models have a larger relative contribution from MTB channel. Since MTB BSSs are typically less massive compared to those created via collisions, the residual 
lifetime for the MTB BSSs is longer than those produced via collisions. Throughout the studied range in p c , our models indicate that a significant fraction of 
today's observed BSSs had been created many Gyrs ago. 



BSS formation channels. In particular we divide all BSSs in 
our models in 5 different categories based on their detailed 
dynamical history. We find that collisions play a major role 
in BSS formation for a large range of central densities typi- 
cal of globular clusters (p c > 10 3 M Q pc -3 ; Figures|4]|5] and 
|6). BSSs in lower density clusters (p c < 10 3 M Q pc~ 3 ) can 
have a significant contribution from the MTB channel (up to 
« 60% in some of our models). This also suggests that for 
open cluster densities (typically < 10 2 M©pc -3 ) the domi- 
nant channel of BSS formation should be MTB. This picture 
can only change if in reality the degree of mixing from the 
MTB channel for some reason is very different than that from 
collisions. One extreme case is, for example, collisions lead 
to minimal mixing and no Hydrogen can be supplied to the 
core. If this is true then there would be no rejuvenation from 
collisions and thus, collisions would produce BSSs that live 
in the BSS stage only a short time. If, however, degree of re- 
juvenation from collisional channels and MTB is similar, then 
our results remain valid. 

The above finding is of great interest. It has been seen in 
the observed BSS populations in the Galaxy that the num- 
ber of BSSs in the core does not s how a positive cor relation 
with the collisional V parameter (iDavies et al.ll2004l) . It has 
also been noted that there is a strong correlation between the 



number of core BS Ss and the binary fraction in the cluster 
( Knig ge~et alj 1200^) . These results had led to conclusions 
in some studies that MTB from undisturbed prim ordial bina- 
ries is the dominant formation channel of BSSs (Solli ma et al.l 
2008). Our results contradict this conclusion and suggest that 
binary mediated collisions is the dominant BSS formation 
channel instead. Although, MTB is not the major contrib- 
utor in BSS formation in typical globular clusters, binaries 
play a crucial role in BSS formation. Because of the higher 
interaction cross-section of the binaries they are much more 
likely to interact compared to the single stars. Moreover, sin- 
gle MS stars tend to reside outside the core at late times due 
to mass-segregation. For typical globular cluster velocity dis- 
persions the branch ing ratio for collisions in binary interac- 
tions is significant (Freg eau et al.l 12003, 2004). BS and BB 
interactions lead to the majority of the collisions in globu- 
lar clusters creating the BSSs (Figure [7] and |8}. Hence, we 
concur with the suggestion that binaries are of great impor- 
tance in creating BSSs in a cluster. However, from our mod- 
els we conclude that this is not because most BSSs are created 
from MTB. The apparent lack of correlation between the spe- 
cific BSS number and the collisional T parameter had been 
noted as suggestive to a predominantly non-collisional chan - 
nel of BSS formation previously (e.g., Sol lima et al.l 2008). 
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Our models suggest that this is clearly not the case. Even for 
the BSSs created via MTB or SE channels, a large fraction of 
them have been dynamically perturbed especially for clusters 
with p c > 10 3 Af Q pc~ 3 . Close inspection of our models re- 
veals that although -/Vbss,c/^Vhb,c do correlate weakly with 
the actual theoretically calculated T in our models, they do 
not with the r o b s estimated taking into account the observa- 
tional difficulties in measuring this quantity. 

We also determined upper limits for the typical t agc for the 
BSSs in our models. In our models the oldest BSSs can be 
almost as old as the cluster itself. The typical ages of the BSSs 
are, however, between 1-5 Gyr (depending on p c ; Figure 
ITOl [TTI and Table [TJ. Thus, present-day BSSs are possibly 
affected by the dynamical history of their host clusters during 
the last several billion years. However, we cautio n that due to 
the assump tion of full-mixing inherent in BSE dHurley et al.1 
l2000l |2002|) . the stellar evolution routine used in our code 
CMC, the i ago values should be tr eated as upper limits, in 
particular for non-rotating stars (e.g., Glebbeek & Pols 2008). 

Throughout this study we have defined BSSs using a sim- 
ple mass-based criterion that is easy in theory but inacces- 
sible in real observations. Although, the BSSs in our mod- 
els identified based on this straight-forward criterion populate 
the expect regions in the HRD, we note that a fraction of the 
theoretically identified BSSs in these models would not have 



been identified as BSSs based on a more observationally mo - 
tivated criteria described, for example, in Leigh et al.l (120071) . 
We have started investigating the properties of the BSSs in our 
models based on this observationally motivated criteria as part 
of another ongoing study. The theoretically extracted BSSs in 
our models is a superset of the BSSs that would be extracted 
using the observationally motivated criterion. 
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Model properties and dynamical properties of BSSs. 
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TABLE 1 

Model properties and dynamical properties of BSSs. 
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NOTE. — t cl denotes the cluster age in Gyr, M denotes the total cluster mass at t cl in 10 4 Mq, f b denotes the overall binary fraction in the cluster, f b c 
denotes the binary fraction in the core, p c denotes the core density in 10 4 Mqpc -3 , and r c denotes the core radius in pc. TVbss an d -^BSS.c denote the 
numbers of BSSs in the whole cluster, and within r c , respectively. The time since formation for the BSSs is denoted by t agc in Gyr. Coll, MTB, and SE denote 
BSSs formed via stellar collisions, stable mass-transfer in a binary, and binary stellar evolution driven mergers, respectively. SS, BS, and BB denote collisions 
via single-single, binary-single, and binary-binary encounters, respectively. Med, and Max denote the median t age and the maximum t agc for each model. 



